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2 s i ABSTRACT 

■ We compare the observed power spectrum of matter found in Papers I and 
Q ■ II with analytical power spectra. We extrapolate the observed power spectra on 

small scales to find the linear power spectrum of matter. We consider spatially 
flat cold and mixed dark matter models with cosmological constant as well as 
open models. We fix the Hubble constant and the baryon density in the middle of 

■ the allowed range and vary the density parameter and the cosmological constant. 
We determine the primordial power spectrum of matter using the power spectrum 
of matter and the transfer functions of analytical models. We take two different 
spectra suggested by observations: one with a sharp maximum at 120/i -1 Mpc and 
a second one with a broader maximum, as found for regions with rich and medium 

& ■ rich superclusters of galaxies, respectively. For both models, the primordial power 

spectra have a break in amplitude; in the case of the spectrum with a sharp 
maximum the break is sharp. We conclude that a scale-free primordial power 
spectrum is excluded if presently available data on the distribution of clusters 
and galaxies represent the true mass distribution of the Universe. 
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1. Introduction 

Here we accept the current inflationary paradigm and assume that the structure evolves 
from Gaussian initial conditions in a Universe dominated by non-baryonic dark matter 
(DM). Under these assumptions the current power spectrum of matter is determined by 
physical processes during the inflationary expansion, and the subsequent radiation- and 
matter-dominated eras. It is believed that the primordial power spectrum, formed during 
the inflation epoch, is a scale-free Harrison- Zeldovich power law. This scale-free power spec- 
trum is transformed during the radiation-dominated expansion, the transformation depends 
on cosmological parameters and on the nature of the DM. Within the presently acceptable 
range of cosmological parameters and possible DM candidates, the transfer function (which 
describes the transformation of the power spectrum) is a smooth function of scale. For this 
reason the present power spectrum should be a smooth function of scale. 

However, recent evidence (Broadhurst et al. (1990), Landy et al. (1996), Einasto et 
al. (1997a), Retzlaff et al. (1998), and Tadros et al. (1998)) indicates that the power spectra 
of galaxies and clusters of galaxies have a spike or peak on scales around I = 120 h~ 1 Mpc 
or wavenumber k = 2ir/l = 0.05 h Mpc" 1 (we designate the Hubble constant as 100 h 
km s" 1 Mpc" 1 ). This scale corresponds to the step size of the supercluster-void network 
(Einasto et al. 1994b, 1997c, hereafter E97c). Atrio-Barandela et al. (1997), and Eisenstein 
et al. (1998) have shown that the peaked power spectrum of galaxies and clusters, and 
the respective angular power spectrum of temperature anisotropy of the cosmic microwave 
background (CMB) radiation are barely compatible with the standard scale-free primordial 
power spectrum in a cold dark matter (CDM) dominated Universe. An agreement is possible 
only for extreme values for cosmological parameters (Hubble constant and baryon fraction 
of the matter) which are almost outside the limits of the range allowed by other data. 

Motivated by the difficulty to reconcile the observed power spectra with CDM-type 
models we try to calculate the primordial power spectrum of matter empirically. We shall use 
the latest determinations of the power spectra for various populations of galaxies and clusters 
of galaxies reduced in amplitude to the power spectrum of matter in the local Universe 
(Einasto et al. 1999a, 1999b, hereafter Papers I and II, respectively). For comparison we use 
theoretical models with cold dark matter (CDM), a mixture of cold and hot dark matter 
(MDM) in spatially flat universe, as well as open models (OCDM). We calculate the transfer 
functions for a set of cosmological parameters chosen in the range of astrophysical interest, 
and find the primordial spectra from theoretical transfer functions and empirical power 
spectra. We calculate for our empirical power spectra the mass function of clusters of galaxies 
using the Press & Schechter (1974) algorithm, and compare them with the empirical cluster 
mass function. 
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2. Comparison with theoretical models 

2.1. Observed power spectra 

We shall use the mean galaxy power spectra determined in Paper I, where we derived 
two mean power spectra based on different observed populations. Clusters of galaxies and 
several galaxy surveys yield a mean power spectrum with a peak on a scale of 120 h~ x Mpc. 
This spectrum, PhdW, characterizes the distribution of clusters and galaxies in a large 
volume that includes rich superclusters of galaxies. For comparison we use also the power 
spectrum found for medium-density regions of the Universe, -Pmd(^) (subscripts HD and 
MD denote high- and medium-density regions). This spectrum is shallower around the 
maximum. The power spectra are determined in the range 0.03 < k < 1 h Mpc -1 . On large 
and intermediate scales (k < 0.2 h Mpc" 1 ) they are found on the basis of 3-dimensional 
galaxy and cluster surveys. Both mean power spectra coincide on small scales, where they 
are based on the reconstruction of the 3-D spectrum from the 2-D distribution of APM 
galaxies. Spectra were reduced to real space and to the amplitude of the power spectrum 
of all galaxies, including dwarf galaxies. This procedure assumes that there exists such a 
mean power spectrum characteristic for all galaxies. Possible errors of this assumption and 
intermediate steps of the data reduction shall be discussed below in Section 3.1. 

Both mean power spectra were reduced in amplitude to the power spectra of matter 
using the bias factor calculated from the amount of matter in voids; the present epoch was 
fixed using the rms density fluctuations in the 8 h~ 1 Mpc sphere, a 8 (Paper II). The biasing 
correction is based on the assumption that the structure evolution in the universe on scales 
of interest is due to gravitational forces alone. Possible errors related to biasing correction 
shall also be discussed below. 

Observed power spectra are shown in Figure 1. They differ only on scales around the 
maximum which is the most interesting part of the spectrum for the present analysis. On 
small scales, k > 0.2 h Mpc" 1 , the power spectra are influenced by non-linear effects. 

2.2. Extrapolation of power spectra on large and small scales 

To compare observed power spectra with theory they are to be reduced to the linear 
case on small scales. To find the linear power spectrum, Peacock (1997) used a method based 
on numerical simulations. We shall use a different method - the comparison with theoretical 
power spectra. Theoretical spectra shall be discussed in more detail in the following Sections. 
Spectra were calculated for CDM and MDM models and for different values of the density 




Fig. 1. — Empirical power spectra of matter, e non-linear power spectrum of matter 

with its 1(7 error corridor (errors were calculated from errors of observed power spectra); 
Phd-Uu-mdm and Pmd-Uti-mdm are the linear power spectra extrapolated on small scales 
using the tilted MDM model with £7 = 0.4, for observed spectra PHD(k) and P MD (k), 
respectively; Puu-cdm is the linear power spectrum extrapolated on small scales using the 
CDM model with Qq = 0.4. On large scales both linear power spectra were extrapolated 
using the standard CDM model with Q = 1- 

parameter, Qq, spectra were COBE normalized. The comparison shows that theoretical 
spectra, calculated for high values of the density parameter, have amplitudes much higher 
on small scales than the amplitude of observed power spectra. The shape of theoretical 
power spectra based on the CDM model is also rather different from the shape of observed 
power spectra on scales between the maximum of the spectrum and the beginning of the 
non-linear part of the spectrum at k pa 0.2 h Mpc -1 . The best fit, both of the amplitude 
(near the beginning of the non-linear scale) and of the shape, is obtained with a tilted MDM 
model with a power index of n = 1.1 on large scales, and a density parameter Qq — 0.4 (see 
Figure 1). For CDM models the best fit is obtained (also of the amplitude and the shape) 
for a model with density parameter Q pa 0.2; but now the fit is good in a much narrower 
scale interval, see Figure 2 below. For comparison we show in Figure 1 an alternative linear 
extrapolation with a CDM model and a density parameter, Q = 0.4. In the following 
discussion we shall use only the extrapolation with the MDM model. Note that for our 
choice of h and VL V (for the MDM model), these parameters corresponds to a neutrino rest 
mass (or the sum of the masses of 2 or 3 stable neutrino species) of pa 3.4 eV. 

Now we consider the extrapolation to large scales. The comparison of observed power 
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spectra with theoretical ones shows that all COBE normalized models with low density 
parameter Q have amplitudes which are higher on large scales than those of the observed 
power spectrum (see Figure 2 below). For this reason we have used the standard CDM 
model with fio = 1 to extrapolate the observed spectra to large scales. This model has the 
lowest amplitude on large scales and yields the best and smoothest extrapolation of observed 
spectra on large scales. On these scales there are no essential differences between theoretical 
spectra of CDM and MDM models. 

We consider these extrapolations on small and large scales as estimates of the true 
linear power spectrum of matter. In the following we shall use the term "empirical power 
spectrum" to denote this linear power spectrum. We recall that it is reduced to real space, 
and its amplitude to the amplitude of the matter density fluctuations. As we have two 
observed mean power spectra of galaxies, characteristic for different populations, we also 
have two empirical power spectra. 

2.3. Theoretical models 

Previous studies (Atrio-Barandela et al. 1997, Eisenstein et al. 1998) have shown that 
it is difficult to make models agree with optical and CMB observations using conventional 
cosmological parameters and a scale-free primordial power spectrum. Here we shall use a 
different approach. We accept cosmological parameters in the middle of the allowed region, 
and try to find possible restrictions to the primordial (initial) power spectrum. 

To calculate the analytical power spectrum we use a Hubble parameter h = 0.6. This 
value is a compromise between the value h = 0.55 ±0.10 suggested by Sandage & Tammann 
(1997) and h = 0.75 ± 0.06 favored by Freedman (1997). We adopt a density of matter 
in baryons of = 0.04 (in units of the critical density). This parameter is in agreement 
with recent nucleosynthesis results fl^h 2 = 0.01 — 0.02 (Songaila, Wampler & Cowie 1997, 
Schramm & Turner 1998, Olive 1997, Buries & Tytler 1998, Steigman 1998, Turner 1998, 
see also Ostriker and Steinhardt 1995). 

We consider models with various density parameter Q Q = Qf> + &dm, as the density 
is presently the parameter which is determined with less accuracy. Local methods of the 
determination of the mass-to-light ratio in clusters and groups of galaxies yield Q 0.2 
(Bahcall 1997). Methods sensitive to the global value of the density parameter yield diverging 
values. Dekel, Burstein & White (1997), using the POTENT method, obtain a value close 
to unity; another treatment of the peculiar velocity field yields a lower value (Freudling et 
al. 1998). The distant supernova project suggests that the Universe is speeding up, i.e. that 
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the vacuum energy contributes about 60 % of the critical density, and that the matter density 
is about 40 % of it, close to results obtained with local methods (Perlmutter et al. 1998, 
Riess et al. 1998). Indirect methods based on the cluster abundance evolution favor a value 
of n = 0.2-0.4 (White, Efstathiou & Frenk 1993, Bahcall, Fan & Cen 1997, Fan, Bahcall & 
Cen 1997, Bahcall & Fan 1998). Further evidence for a low-density Universe comes from the 
analysis of the Lyman-alpha forest by Weinberg et al. (1998), who found for a flat Universe 
fi = 0.34to;o9 (1°" errors). A recent analysis of all available data by Turner (1998) favors 
h = 0.60 - 0.65, Q « 0.4 and Q A « 0.6. 

To fix the amplitude of the analytical power spectrum on large scales we use the four- 
year COBE normalization (Bunn & White 1997). We cannot use simultaneously the a 8 
normalization on small scales as we do not know which model best fits both normalizations. 
It is just our goal to try to find a model which satisfies both normalizations. 

We use three sets of models: spatially flat models with cold and mixed dark matter, and 
open models with cold dark matter. In first two models we use the cosmo logical constant. Its 
necessity is strongly favored by recent data on cluster abundance and supernova explosions 
at high redshifts (Bahcall et al. 1997, Perlmutter et al. 1998, Riess et al. 1998). We vary 
the density parameter between Q = 0.2 and Qo — 1, and we choose the cold dark matter 
(CDM) content to obtain a spatially flat model, f2& + + = 1- Models are calculated 
for a number of spectral indices on large scales, n = 1.0, 1.1, .. . 1.4. 

We use a similar set of parameters for mixed dark matter (MDM) models with the 
only difference that a hot dark matter component with density parameter Vt v = 0.1 was 
added. The density parameter of the cold dark matter was decreased to get a spatially flat 
model. We use a small fraction of the hot dark matter for two reasons. First, most direct 
and indirect methods suggest that f2 should be less than unity. Q v should be chosen in 
agreement with the relation Q v /Qq ~ 0.2 required to get the correct value of as (see, e.g., 
Pogosyan & Starobinsky (1993), and Pogosyan & Starobinsky (1995) for the case of a tilted 
initial spectrum with n > 1). Second, in order to build up the fine filamentary structure of 
faint galaxies observed in supervoids (Lindner et al. 1995), the fraction of matter in the cold 
dark matter must be significantly larger than that of the hot component (Frisch et al. 1995). 

In the open models (OCDM) we used values of f2 = 1.0, 0.9, ... 0.3 for the matter 
density. The baryon density was fixed as in all other models, and the density of the CDM 
was chosen appropriately. 

Theoretical power spectra for our models are plotted in the left panels of Figure 2 
together with the empirical linear power spectra of matter. Models were calculated with 
the CMBFAST package of Seljak & Zaldarriaga (1996). We shall discuss the linear power 
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spectra and the comparison of model spectra with observations in subsequent sections below. 

2.4. Primordial power spectra 

The CMBFAST package yields the transfer function for each of our models. We have 
used this function to calculate primordial power spectra, 

P mit (k)=P(k)/T 2 (k), (1) 

where T(k) is the transfer function. In the right-hand panels of Figure 2 we plot the ratio of 
the primordial power spectrum to the scale-free primordial power spectrum, Pmit(k) / Po(k) , 
where Po(k) ~ k. We plot here only results found for the peaked empirical power spectrum 
-Pffn(^)- A comparison for the flatter power spectrum P M£) (fc) (which represents medium 
rich regions) is given in Figure 3. In both cases we have used the empirical power spectra 
extrapolated on small and large scales as described above. The discussion of primordial 
power spectra follows in Section 3.3. 

3. Discussion 

3.1. Power spectra 

We compare now empirical power spectra with theoretical ones on large, intermediate 
and small scales. 

On large scales the amplitude of COBE normalized theoretical power spectra depends on 
the density parameter. Low-density CDM and MDM models dominated by the cosmological 
constant have amplitudes higher than models of critical density. Figure 2 demonstrates 
that on large scales empirical power spectra have amplitudes which are compatible with 
amplitudes of COBE normalized theoretical power spectra only for high-density CDM and 
MDM models with Jl ~ 1. For this reason a smooth extrapolation of empirical power 
spectra on large scales was possible only for high-density models. The highest amplitude 
on large scales (for extrapolation of observed power spectra) which is still in satisfactory 
agreement with the COBE normalization of theoretical spectra and with the upper limit of 
the error corridor of empirical spectra, is the CDM or MDM spectrum with Qq = 0.7 (the 
CDM and MDM model spectra are similar in this range of scales). 

On intermediate and small scales, k > 0.05 h Mpc -1 , the amplitude of theoretical 
power spectra of CDM and MDM models varies considerably with the density parameter. 
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Only theoretical power spectra with a low density parameter have amplitudes similar to the 
amplitude of the empirical power spectra. For CDM model the best agreement is achieved for 
a model with Q « 0.2, and for an MDM model for Q 0.4 (see Figure 2). On these scales 
CDM and MDM models are in agreement with other evidence, which suggests that high- 
density models with Qq w 1 are excluded by a large margin (Ostriker and Steinhardt 1995, 
Bahcall, Fan, Cen 1997, Fan, Bahcall, Cen 1997, Borgani et al. 1997, Bahcall & Fan 1998, 
Turner 1998, Weinberg et al. 1998). Low-density theoretical power spectra have, however, 
maxima at k < 0.01 h Mpc -1 with an amplitude of the spectrum which is much higher than 
that of the observed peak at k — 0.05 h Mpc -1 . Available data do not support the presence 
of a rising spectrum on these scales: spectra found from cluster and galaxy data decrease in 
amplitude toward large scales. 

Open models have a smaller increase of the amplitude on large scales for a wide range of 
the density parameter, and are thus in a better agreement with a reasonable extrapolation 
of the observed power spectrum on large scales. This agreement is lost around the maximum 
of the empirical power spectrum, and both versions of empirical power spectra deviate here 
from OCDM models by a large margin. 

This comparison shows that it is impossible to satisfy the shape of the empirical power 
spectrum with models with a fixed density parameter simultaneously on large and small 
scales. Models which fit empirical spectra on large scales are incompatible with empirical 
spectra on small scales and vice versa. This is the main conclusion obtained from the 
comparison of cosmological models with the data. 

3.2. Primordial power spectra 

Now we discuss primordial spectra derived from the empirical power spectrum, Phd{^)- 
On large scales there is little difference between CDM and MDM models. The main feature of 
primordial power spectra is the presence of a spike at the same scale as that of the maximum 
of the empirical power spectrum. On scales shorter than that of the spike, the primordial 
spectrum can be well approximated by a power law (a straight line in a log-log plot). The 
slope of this line varies with the cosmological density parameter accepted for theoretical 
models. For models with high cosmological density, f2 > 0.5, the slope is negative, n < 1; 
for models with low density, Qo < 0.5, it is positive, n > 1. Such tilted primordial power 
spectra (i.e. spectra with n ^ 1) have been used to increase the accuracy of the model. 
This approximation breaks down, if we consider the whole scale interval. On large scales the 
primordial power spectrum can also be approximated by tilted models. The power index of 
the approximation is, however, completely different from the index suitable on small scales, 
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as seen in Figure 2. 

Additionally, there is a considerable difference in the amplitude of the primordial power 
spectrum on small and large scales. For most values of the cosmological density parameter 
the amplitude on small scales is lower than on large scales (compared to the scale-free 
primordial spectrum). For low values of fio the effect has the opposite sign: the amplitude 
of the primordial spectrum on small scales is higher than on large scales. Such an effect 
has been noticed already by Lesgourgues et al. (1998). The amplitude of the break depends 
on the model used, and in MDM models it is larger than in CDM models. The primordial 
power spectrum for open model has similar features, only its amplitude on large scales rapidly 
increases if we consider more open models. This is due to the fact that the amplitude of 
theoretical power spectra of these models on scales k > 0.05 h Mpc -1 is much lower than 
the amplitude of the empirical power spectrum. 

Primordial power spectra, which are found from the comparison of the flat empirical 
power spectrum, Pml>(&), with theoretical spectra, are plotted in Figure 3. On very large and 
small scales both empirical power spectra are identical, and consequently primordial power 
spectra are identical, too. Differences exist in the medium scale range. Here primordial 
power spectra have a smoother transition from long to small scales with no sharp spike 
at wavenumber k = 0.05 h Mpc -1 . However, a change of the spectral index and amplitude 
around this wavenumber (a "break" ) is very well seen. This property of the primordial power 
spectrum is similar for spectra derived on the basis of CDM, MDM and OCDM models. An 
independent analysis of the initial power spectrum was made by Adams, Ross & Sarkar 
(1997) using the APM galaxy power spectrum (which is identical to our Pmd(&))- The 
presentation of the initial power spectrum (their Fig. 1) is slightly different from ours, the 
main features of the spectrum are, however, well seen, with similar conclusions. 

Previous studies have shown that it is practically impossible to bring the classical models 
into agreement with new data on the power spectrum by varying cosmological parameters. 
Our analysis confirms these results. The reason for the discrepancy lies in the shape of the 
empirical power spectrum: it is much narrower than theoretical spectra. A change of the bias 
parameter or a normalization on large scales does not change these conclusions: they would 
only shift the features present in Figures 2 and 3 up or down, without changing considerably 
the shape of the primordial spectrum. 
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3.3. Models with low density parameter 

The previous analysis has shown that on scales k > 0.05 h Mpc -1 the empirical power 
spectrum can be well approximated by a MDM model with density parameter Q Q fa 0.4. A 
slightly less accurate approximation in a narrower scale interval is provided by a CDM model 
with Q ~ 0.2. An open model with density parameter Q fa 0.4 gives an approximation 
which fits the empirical power spectrum Pmd both on large and small scales (but not on 
intermediate scales). A spike or break of the primordial spectrum at k — 0.05 h Mpc -1 could 
be avoided using one of these models, if we allow for the higher amplitudes on large scales 
of the empirical power spectrum (for CDM and MDM models), or a shallower turnover on 
medium scales (for the OCDM model). Examples of such spectra for MDM and OCDM 
models are shown in Figure 4. It is reasonable to investigate the possibilities of such models. 

There are two major reasons to reject these model power spectra. First, the directly 
observed power spectrum does not have a rising section above the maximum at k = k = 
0.05 h Mpc -1 , needed to fit low-density CDM or MDM models, or a very shallow maximum 
needed for the low-density OCDM model. Second, as shown by numerical experiments 
discussed in Paper I, models with rising or shallow power spectra in this scale interval have a 
distribution of superclusters which is very different from the distribution of real superclusters. 
In a model with a broad power spectrum, such as all scale-free CDM, MDM and OCDM 
models, the distribution of rich superclusters is close to a random one. Real rich superclusters 
form a quasi-regular lattice, which surrounds large voids of diameter fa 120 h^ 1 Mpc (E97c). 
A quantitative measure of the regularity is the cluster correlation function on large scales, 
which oscillates with a period that corresponds to the step size of the lattice (Einasto et 
al. 1997b, E97b). 

To test model power spectra for the regularity of superclusters we have calculated via 
Fourier transformation correlation functions for models shown in Figure 4. To compare with 
the correlation function of clusters of galaxies in rich superclusters we multiplied correlation 
functions with factors which correspond to the mean difference in the amplitude of power 
spectra of matter and clusters in rich superclusters. Correlation functions are plotted in 
the right panel of Figure 4. We see that the correlation function calculated from the power 
spectrum P HD is oscillating with a low amplitude, both MDM and OCDM models do not 
show signs of oscillations. We note that the amplitude of oscillations of the observed cluster 
correlation function is exaggerated due to the elongated form of the observed cluster sample. 
However, as shown by simulations, it is impossible to generate regular oscillations of the 
cluster correlation function from a random distribution of superclusters. 

To estimate the statistical significance of the difference between observed and model 
power spectra we can apply the procedure used by E97b. The broad power spectrum of MDM 
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and OCDM models can be approximated by a geometrical model with randomly distributed 
superclusters. E97b generated 1000 simulated catalogs of clusters, applied selection criteria 
similar to those used in real cluster samples, selected very rich superclusters, and calculated 
the respective correlation functions. Very rich superclusters cause peaks and valleys in the 
correlation function; the oscillatory behavior of the correlation function was characterized 
by the period and amplitude of oscillations, their rms scatter, and by a parameter which 
characterizes the regularity of oscillations. This test shows that some individual combinations 
of oscillation parameters are in agreement with respective parameters of the real cluster 
correlation function with a probability of about 1 %, but a simultaneous occurrence of all 
parameters has a much smaller probability (of the order of 10" 6 ). This test shows that there 
is a statistically highly significant difference between samples which generate broad and 
peaked power spectra. A similar conclusion has been obtained by Luo & Vishniac (1993) 
by the comparison of the 1-D power spectrum found for the Broadhurst et al. (1990) pencil 
beam survey with the theoretical spectrum calculated on the basis of the standard theory. 
Significant differences between cluster power spectra and spectra based on various variants 
of CDM models have been noticed by Retzlaff et al. (1998) and Tadros et al. (1998); similar 
differences have been found also in case of the power spectrum based on the 2-D distribution 
of APM galaxies by Peacock (1997) and Gaztahaga & Baugh (1998). 



3.4. Theoretical implications 

A consequence of the broken-scale-invariant (BSI) primordial power spectrum is that 
one cannot draw conclusions based on the exact scale-free primordial spectrum, such as 
the use of the zero crossing of the correlation function, to determine the density parameter 
through the parameter T = Qoh. Another conclusion is that one should not expect the slopes 
n,s and nx of power spectra of primordial matter and gravitational waves to be independent 
on k, especially, around the critical point k = 0.05 h Mpc -1 . Here ns and ut correspond to 
matter (scalar perturbations) and gravitational waves (tensor perturbations), respectively. 
In the latter region, differential relations between these quantities following the slow-roll 
approximation during inflation may become invalid. 

Also, we have to understand the reason why the primordial power spectrum deviates 
from the classical Harrison- Zeldovich spectrum. A very general way to obtain such a behavior 
requires some kind of phase transition which occured during inflation, about 60 e-folds before 
the end of the inflation era. There is one specific model of this type which has an exact 
solution. It was suggested by Starobinsky (1992), and leads to a definite BSI spectrum. 
It was recently confronted with observational data by Lesgourgues et al. (1998). However, 
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there exist a number of other possibilities (see recent discussions by Adams, Ross & Sarkar 
(1997) and Starobinsky (1998)) which should be explored, too. 

4. Conclusions 

The main goal of this series of papers was to determine the mean empirical power 
spectrum of galaxies, to reduce it to the matter, and to compare with theoretical models 
to see whether it is possible to fit empirical spectrum with models. This procedure rests 
on several assumptions and conclusions based on these assumptions. The main assumptions 
are: 

(1) there exists a power spectrum which represents all galaxies in a fair sample of the 
Universe (including faint dwarf galaxies); 

(2) the dynamical evolution of the Universe on scales of interest is determined by gravity 

only; 

(3) density perturbations grow from small random fluctuation generated in the early 
stage of the evolution; 

(4) the dynamics of the Universe is dominated by cold dark matter with some possible 
mixture of hot dark matter; 

(5) galaxy samples of various environment, morphology and luminosity can be approx- 
imated by particles in numerical simulations chosen in certain threshold density intervals. 

The first assumption is based on observation of power spectra of galaxy samples with 
different absolute magnitude limits: spectra are identical in shape and amplitude if the 
sample contains sufficiently faint galaxies (GE92). The scatter of observed data points of 
power spectra of faint galaxy samples is about 10 % which can be attributed to the cosmic 
scatter. Thus we can expect that the assumption itself does not introduce any systematic 
error to the analysis. The last assumption is based on our experience in comparing real and 
simulated galaxies, collected in last ten years; its justification was analyzed in detail in Paper 
II of this series. 

The second, third and fourth assumptions form main paradigms of modern cosmol- 
ogy. The second and third assumptions are critical in the whole chain of reduction steps 
of observed power spectra. Under these assumptions, due to gravitational instability, the 
structure evolution in high- and low-density regions is different: in low-density regions mat- 
ter remains in the primordial form, whereas in high-density regions matter collapses and 
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forms galaxies and systems of galaxies. A clump of matter collapses if its density is high 
enough, thus the gravitational character of the evolution leads to the conclusion that galaxy 
formation is a threshold phenomenon - there exists a critical threshold density, needed for a 
clump of matter to collapse. Here the actual density of the primordial matter is important. 
Groups and clusters of galaxies have a characteristic scale of the order of 1 h~ x Mpc, thus 
calculated densities must not be smoothed over scales exceeding this value. Various galaxy 
populations can be simulated using different threshold densities (5th assumption). A thresh- 
old density, equal to the mean density of matter, divides the population of all galaxies from 
that of matter in voids, higher threshold densities correspond to luminous galaxies, and to 
galaxies in clusters. 

The division of matter into unclustered and clustered populations affects the respective 
power spectra in a simple way, so that power spectra of galaxies, systems of galaxies, and 
of matter are similar in shape. We can reduce spectra of different populations to the power 
spectrum of all galaxies by a shift in the amplitude only. Simulation show that relative errors 
of this reduction procedure are of the order of 1 % for galaxies and about 5 % for clusters. 

The comparison of power spectra of simulated galaxy and cluster samples in real and 
redshift space has shown that the influence of peculiar velocities in groups and clusters is 
serious on small scales only. To avoid complications with the reduction of galaxy power 
spectra to real space, we have accepted on small scales (k > 0.1 h Mpc -1 ) the power 
spectrum of APM galaxies as the power spectrum of all galaxies in real space. This power 
spectrum is deduced from 2-D distribution of galaxies, and has no redshift distortions. Here 
we assume that the APM sample is a fair sample of all galaxies in the Universe on these 
scales. The main error of the amplitude of the power spectrum of all galaxies is related with 
this assumption. The amplitude may have a relative error of the order of 15 %, which is the 
largest possible error of the power spectrum of all galaxies. 

On large scales the mean power spectrum is determined essentially by clusters of galax- 
ies. Numerical simulation confirm that on these scales the redshift correction (due to the 
contraction of superclusters) affects only the amplitude of power spectra; possible systematic 
errors involved are negligible in comparison with sampling errors. 

The final step in deriving the empirical power spectrum of matter is the reduction of the 
mean power spectrum of all galaxies to that of matter (the conventional biasing correction). 
Here we use again our basic assumptions that the structure evolution is due to gravity 
alone, and that initial density fluctuation have a Gaussian distribution. Additionally we use 
our last assumption that galaxy populations can be approximated by samples of simulated 
particles chosen in various threshold density intervals. As discussed above, under these 
assumptions the matter power spectrum can be found by shifting the galaxy spectrum in 
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amplitude. The shift is determined by the fraction of matter in the clustered population. 
There are two sources of error in this procedure: the uncertainty of the threshold density 
level, Qo, which divides the low-density population of primordial particles in voids from the 
clustered population of particles associated with galaxies; and errors due to uncertainty of 
cosmological parameters which define the speed of void evacuation and the fraction of matter 
in high-density regions (see Paper II for details). The relative error of the threshold density 
level was estimated by Einasto et al. (1994a) to be 10 %, which leads to a 3 % error of the 
fraction of matter in the clustered population and respective biasing parameter. The second 
error can be estimated from data given in Paper II (see Table 2 and Figure 4), it is 7 %. 
Additional sources of error are the fuzziness of the threshold density, go, and the assumption 
that the distribution of luminous galaxies in superclusters, clusters and groups is similar to 
the distribution of all matter. The analysis done in Paper II shows that these errors do not 
influence the shape of the power spectrum in the scale range of interest (corresponding errors 
are of the order of 1 %); the errors in the amplitude are also small, less than or equal to 5 %. 
The total rms relative error of the biasing factor is 10 % (all errors are la errors). 

To summarize the overview of the data reduction procedure we can say that the gravita- 
tional character of the structure evolution implies that intermediate steps of data reduction 
introduce no noticeable error to the shape of the power spectrum, possible errors are in the 
amplitude of the power spectrum only. The total rms error of the amplitude is essentially 
determined by two errors: the error of using the amplitude of the APM galaxy power spec- 
trum as the amplitude of the power spectrum of all galaxies in a fair sample of the Universe; 
and the error due to the uncertainty of model parameters used in the reduction of the galaxy 
power spectrum to matter. 

Possible errors of the shape of the power spectrum are intrinsic, due to differences of 
power spectra of different populations. To quantify possible differences in the shape of the 
power spectrum of matter we have formed two mean power spectra, characteristic for samples 
which cover regions of the Universe including high-density superclusters, and medium-density 
superclusters only. Differences between power spectra of these samples may partly be due 
to unknown errors of data reduction used in various samples, or to the geometry of samples 
(thin slices in case of the LCRS). 

Empirical mean power spectra of matter are extrapolated on small scales (to get linear 
power spectra) and on large scales (to have the spectrum defined on all scales for comparison 
with theory). Power spectra found for high- and medium-density regions are identical on 
small scales (k > 0.1 h Mpc -1 ), and on large scales (k < 0.01 h Mpc^ 1 ). On medium 
scales they are different. The power spectrum found for regions which include high-density 
systems with rich superclusters has a peak at k — 0.05 h Mpc^ 1 , and an almost constant 
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power index n = —1.9 for scales in the range 0.05 < k < 0.2 h Mpc -1 . The power spectrum 
of medium-density regions has a shallower shape around the maximum, the amplitude near 
the maximum is lower by a factor of 2. 

Empirical power spectra can be well approximated by theoretical model spectra based 
on various dark matter models in a limited scale range. The best approximation for the 
empirical power spectrum for high-density regions on medium and small scales is provided by 
a spatially flat MDM model with cosmological density parameter Q = 0.4, baryonic density 
Qb = 0.04, cold dark matter density VLcdm = 0.26, hot dark matter density Q nu = 0.1, 
and cosmological constant term Q\ = 0.6; Hubble parameter was fixed at h = 0.6. This 
approximation breaks down on large scales: the model power spectrum continues to rise on 
scales k < 0.05 h Mpc -1 , whereas the amplitude of the empirical power spectrum decreases. 
The impossibility to represent empirical power spectra with one theoretical model is one of 
the main results of this series of papers. 

If we make use of our fourth assumption that the dynamical evolution of the Universe 
is dominated by some sort of dark matter, then we can use our empirical power spectra to 
calculate the primordial power spectrum. These calculations show that for both variants 
of the empirical power spectra the primordial power spectrum has a break, i.e. it can be 
approximated on large and small scales by two power laws with different power indices. This 
is due to the fact that it is impossible to approximate the empirical power spectrum with 
one single model. The strength of the break depends on the cosmological parameters used 
for theoretical models, and the shape of the break is different for empirical power spectra 
found for high- and medium-density regions. 

Our conclusions are based on the tacit assumption that galaxy and cluster samples 
used represent the true matter distribution in the Universe. To check this assumption new 
deeper redshift data are needed. The best dataset to make this crucial check is the red- 
shift survey of giant elliptical galaxies planned in the Sloan Digital Sky Survey (SDSS, see 
astro-ph/ 98101301 , |astro-ph/980"9085| , |astro-ph/9809T79| , |astro-ph/980"53ll ). Giant elliptical 



galaxies are concentrated to clusters and superclusters, and their distribution on large scales 
provides the best test for the shape of the power spectrum of matter on scales around the 
maximum. Presently we can say that the possibility of a broken-scale-invariant primordial 
power spectrum deserves serious attention. 
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Appendix: Cluster mass function 

We use the empirical power spectra of mass to calculate the cluster mass function. The 
latter can be compared with the observed cluster mass function which allows us to check 
both sets of data for consistency. 

We calculated the cluster mass function using the algorithm, developed by Suhhonenko 
& Gramann (1998) on the basis of the Press-Schechter (1974) theory. This function provides 
an important constraint for the power spectrum (White, Efstathiou & Frenk 1993, Cen 
1998). The cluster mass function was found for both empirical power spectra, -Phx>(/c), and 
-Pmd(^), reduced to the linear regime. Mass functions were derived for three values of the 
density parameter, Q = 1.0, 0.6, 0.4. Results are shown in Figure 5, together with the 
observed cluster mass function according to Bahcall & Cen (1993). Our calculations show 
that both variants of the empirical power spectra yield a cluster mass function which has 
a shape close to the shape of the observed cluster mass function. The differences between 
mass functions for spectra Phd^) an d -Pmd(^) are small, in other words, the mass function 
is not sensitive to such differences of power spectra. This is not surprising since clusters of 
galaxies are formed by density perturbations on 1 - 10 h' 1 Mpc scales where both empirical 
power spectra are identical. On the other hand, the amplitude of the mass function is very 
sensitive to the density parameter of the Universe used for calculations. The amplitude 
is evidently too high if one assumes a high density parameter of the Universe. The best 
agreement is obtained for the case of Q = 0.6. A value of Q = 0.4, as favored by cluster 
abundance evolution constraints (Bahcall & Cen 1993, Bahcall, Fan & Cen 1997), yields a 
lower amplitude of the calculated cluster mass function. 

The focus of this series of papers is on the empirical power spectrum and on its conse- 
quences for the primordial power spectrum. We do not consider the calculation of the cluster 
mass function for empirical power spectra as a density determination. The only conclusion 
we can draw from this comparison is that our power spectra yield cluster mass functions in 
good agreement with the observed cluster mass function. A relatively high value of the den- 
sity parameter favored by this test can be explained if we assume that the biasing correction 
found from the fraction of matter in high-density regions is too large (see Paper II). 
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Fig. 2. — The empirical linear power spectra of matter compared with theoretical and pri- 
mordial power spectra. Left: present power spectra; right: primordial spectra; in upper 
panels for CDM models, in middle panels for MDM models, and in lower panels for OCDM 
models. Solid bold lines show the empirical linear matter power spectrum found for regions 
including rich superclusters, Phd^); dashed bold lines show the empirical linear power 
spectrum of matter Pmd(^) for medium dense regions in the Universe. Model spectra with 
Qq = 1.0, 0.9, . . . 0.2 are plotted with thin lines; for clarity the models with Q = 1.0 and 
f2o — 0.5 are drawn with dashed lines. In the case of MDM the model of lowest density 
parameter is that with Qq = 0.25. Only spectra with spectral index n = 1 on large scales 
are plotted. Primordial power spectra are shown for the peaked matter power spectrum, 
-Pffo(^); they are divided by the scale-free spectrum, P(k) ~ k. 



- 21 - 



MDJniM.OO 
MD_init_0.50 




10 



- P_MDJnit_1.00 
P_MD_init_0.50 



D_ 




0.1 
k [h Mpc"'] 



0.1 
k [h Mpc"'] 



Fig. 3. — Primordial power spectra for the empirical linear matter power spectrum, Pmd(^); 
divided by the scale-free spectrum, P(k) ~ k. Left: CDM model with decreasing density 
parameter Q from below on the right side of the panel; right: MDM with decreasing Q 
from below on the right side of the panel. Models with Q = 1.0 and Qo = 0.5 are plotted 
with dashed lines for clarity. 
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Fig. 4. — Power spectra (left) and correlation functions (right) of MDM and OCDM mod- 
els with density parameter Q = 0.4, compared with the linear empirical power spectrum 
of matter and correlation function of clusters of galaxies in rich superclusters. The power 
spectrum of the MDM model is calculated with spectral index n = 1.1. Cluster correla- 
tion functions are calculated via Fourier transform from power spectra of matter, and are 
enhanced in amplitude by a biasing factor 7.7 which corresponds to the mean difference 
between respective power spectra. The observed cluster correlation function is the one for 
Abell-ACO clusters in very rich superclusters as derived by E97b. 
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Fig. 5. — Cluster mass functions. Open and filled circles with error bars indicate the observed 
cluster mass function according to Bahcall and Cen (1993); solid and dashed lines show mass 
functions calculated for the linear empirical power spectra, -Phx>(/c) and Pmd(^), for density 
parameters f2 = 1-0, 0.6, 0.4. 



